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A measurement is presented of differential cross sections for Higgs boson (H) produc¬ 
tion in pp collisions at ^/s = 8 TeV. The analysis exploits the H —> 77 decay in data 
corresponding to an integrated luminosity of 19.7 fb^^ collected by the CMS experi¬ 
ment at the LHC. The cross section is measured as a function of the kinematic prop¬ 
erties of the diphoton system and of the associated jets. Results corrected for detector 
effects are compared with predictions at next-to-leadrng order and next-to-next-to- 
leading order in perturbative quantum chromodynamics, as well as with predictions 
beyond the standard model. For isolated photons with pseudorapidities \r]\ < 2.5, 
and with the photon of largest and next-to-largest transverse momentum (pj) divided 
by the diphoton mass myj satisfying the respective conditions of pj/niyy >1/3 and 
>1/4, the total fiducial cross section is 32 ± 10 fb. 
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1 Introduction 

In 2012, the ATLAS and CMS Collaborations announced the observation ||l]-[3l of a new boson 
with a mass of about 125 GeV, with properties consistent with expectations for the standard 
model (SM) Higgs boson. The Higgs boson (H) is the particle predicted to exist as a conse¬ 
quence of the spontaneous symmetry breaking mechanism acting in the electroweak sector of 
the SM This mechanism was suggested more than fifty years ago, and introduces a com¬ 
plex scalar field, which gives masses to W and Z bosons HTj- ITTI . The scalar field also gives 
mass to the fundamental fermions through a Yukawa interaction |j5l. Couplings and spin of 
the new boson are found to be consistent with the SM predictions lfT2] - [l^ . A measurement of 
the pp —t H —> 77 differential cross section as a function of kinematic observables investigates 
possible deviations in distributions related to production, decay, and additional jet activity. It 
provides a check of perturbative calculations in quantum chromodynamics (QCD), and can 
point to alternative models in the Higgs sector. A similar analysis has been carried out by the 
ATLAS Collaboration in diphoton and four-lepton decay channels IITtI - IT^ . 

Despite its small branching fraction of ^0.2% Il20]| predicted by the SM, the H —> 77 decay 
channel provides a clean final-state topology and a precise reconstruction of the diphoton mass. 
The dominant background arises from irreducible direct-diphoton production and from the re¬ 
ducible pp —> 7 +jets and pp ^ jets final states. The relatively high efficiency of the H ^ 77 
selection makes this final state one of the most important channels for observing and investi¬ 
gating the properties of the new boson. 

In this paper, the cross section is measured as a function of the kinematic properties of the 
diphoton system, and, in events with at least one or two accompanying jets, also as a function 
of jet-related observables. Two isolated photons are required to be within pseudorapidities 
I?/1 < 2.5, and the photon with largest and next-to-largest transverse momentum (pj) must 
satisfy the respective conditions of pj/nijj >1/3 and >1/4, where niyy represents the dipho¬ 
ton mass. The transverse momentum p^^ and the rapidity of the Higgs boson, observ¬ 
ables related to the opening angle between the two photons, and the number of jets Njets with 
Pt > 25 GeV produced in association with the diphoton system, are defined in this inclusive 
fiducial selection. A departure relative to the SM-predicted angular distributions would be an 
important observation, as it could reflect different spin and parity properties 1211 than expected 
in the SM. 

The variables defined with at least one accompanying jet are sensitive to the transverse Lorentz 
boost of the diphoton system. A modification in the corresponding distributions or in the pj^ 
spectrum could signify new contributions to gluon-gluon fusion production of the Higgs boson 
(ggH) 1 22 1 . The variables defined by requiring at least two accompanying jets are related to 
production of H —> 77 through vector boson fusion (VBF); however, given the low event yield 
after selecting two jets, no other selection is applied to enhance this production mechanism. 
This is different from what was done in Ref. liT2l , where an attempt was made to classify the 
events according to the production mechanism. The differential cross sections are therefore 
mainly sensitive to the dominant ggH production mode of the Higgs boson. 

The data correspond to an integrated luminosity of 19.7fb^^ collected at the GERN LHG by 
the GMS experiment in proton-proton collisions at = 8 TeV. The trigger requirements and 
vertex determination are identical to those of Ref. fl2l , while photon selection and event classi¬ 
fication are modified to reduce their dependence on pj and providing thereby a less model- 
dependent measurement. Photons are identified using a multivariate classifier that combines 
information on distributions of shower and isolation variables designed to be independent 
of pj and p'^. The signal yield is extracted by fitting the distribution simultaneously in 
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2 The CMS detector 


all bins of the observables. To improve the sensitivity of the analysis, the selected events are 
categorized using an estimator of the mass resolution that is not correlated with niryy, which 
simplifies the description of the background. Measured distributions are unfolded for detector 
effects and compared to distributions at the generator level from the latest Monte Carlo (MC) 
predictions. 

The paper is organized as follows. After a brief description of the CMS detector and event re¬ 
construction given in Section]^ and of the simulated samples in Section]^ the photon selection 
and event classification are detailed in Section]^ where we also describe the kinematic observ¬ 
ables. Section|^provides the statistical methodology for extracting the signal, and gives details 
on modelling signal and background, and on the unfolding procedure. Systematic uncertain¬ 
ties are detailed in Section Unfolded results are then compared with theoretical predictions 
in Section]^ and a brief summary is given in Section]^ 

2 The CMS detector 

A full description of the CMS detector, together with a definition of the coordinate system and 
the relevant kinematic variables, can be found in Ref. Il23ll . Its central feature is a supercon¬ 
ducting solenoid, 13 m in length and 6 m in diameter, which provides an axial magnetic field of 
3.8 T. The core of the solenoid is instrumented with trackers and calorimeters. The steel flux- 
return yoke outside the solenoid is equipped with gas-ionisation detectors used to reconstruct 
and identify muons. Charged-particle trajectories are measured using silicon pixel and strip 
trackers, within \r]\ < 2.5. A lead tungstate crystal electromagnetic calorimeter (ECAL) and 
a brass and scintillator hadron calorimeter (HCAL) surround the tracking volume and cover 
the region \r]\ < 3. The ECAL barrel extends to \rj\ < 1.48, while the ECAL endcaps cover 
the region 1.48 < \ri\ < 3.0. A lead and silicon-strip preshower detector is located in front 
of each ECAL endcap in the region 1.65 < | 7 /| < 2.6. The preshower detector includes two 
planes of silicon sensors that measure the transverse coordinates of the impinging particles. A 
steel and quartz-fibre Cherenkov calorimeter extends the coverage to |?/| < 5.0. In the (?/,</>) 
plane, for |?/| < 1.48, the HCAL cells map onto 5x5 ECAL crystal arrays to form calorimeter 
towers projecting radially outwards from points slightly offset from the nominal interaction 
point. In the endcap, the ECAL arrays matching the HCAL cells contain fewer crystals. To 
optimize the energy resolution, the calorimeter signals are calibrated and corrected for several 
detector effects |24|. Calibration of the ECAL uses the (^-symmetry of the energy flow, photons 
from 71° —>• 77 and rj —?■ 77 , and electrons from W —> ev and Z —> e^e^ decays. Changes in 
the transparency of the ECAL crystals due to irradiation during the LHC running periods and 
their subsequent recovery are monitored continuously and corrected, using light injected from 
a laser system. 

Photons are reconstructed from clusters of energy deposition into so-called "superclusters" in 
the ECAL ||25| . Events are selected using triggers requiring two photons with different thresh¬ 
olds in energy in the transverse plane (Ej)^ respectively, Ej > 26 and > 18 GeV for the leading 
and subleading photons, and through other complementary selections. One selection requires 
a loose calorimetric identification based on the distribution in energy in the electromagnetic 
cluster, and loose isolation requirements on photon candidates. The other selection requires 
a photon candidate to have a high value of Rg variable, defined as the sum of the energies 
deposited in the array of 3 x 3 crystals centred on the crystal with highest energy deposition 
in the supercluster, divided by the energy of the supercluster. Photons that convert to e+e^ 
pairs before reaching the calorimeter tend to have wider showers and smaller values of Rg 
than unconverted photons. High trigger efficiency is maintained by having both photons sat- 
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isfy either selection. The measured trigger efficiency is 99.4% for events satisfying the diphoton 
preselections described in Section]^ 


3 Monte Carlo samples 

The MC simulation of detector response employs a detailed description of the CMS detector, 
and uses Geant4 version 9.4.p03 Simulated events include additional pp collisions that 
take place in or close to the time span of the bunch crossing, and overlap the interaction of inter¬ 
est. The probability distribution of these pileup events is weighted to reproduce the observed 
number of interactions in data. 

The MC samples for ggH and VBF processes use the next-to-leading order (NLO) matrix el¬ 
ement generator POWHEG (version 1.0) Il27l - |3T|| interfaced with PYTHIA 6.426 1321 . The CTIO 
1331 set of parton distribution functions (PDF) is used in the calculation. The POWHEG gener¬ 
ator is tuned following the recommendations of Ref. Il34l and reproduces the Higgs boson px 
spectrum predicted by the HqT calculation Il35ll3^ . The PYTHIA 6 tune Z2* is used to simulate 
the hadronization and underlying event in pp collisions at 8 TeV. The Z2* tune is derived from 
the Z1 tune 1371 , which uses the CTEQ5L PDF, whereas Z2* adopts CTEQ6L1 I38l . The cross 
section for the ggH process is reduced by 2.5% for all values of nin to accommodate its interfer¬ 
ence with nonresonant diphoton production 11391 . The PYTHIA 6 generator is used alone for the 
VH (where V represents either the W or Z boson) and ttH processes with the CTEQ6L1 PDF 
138 1 and Z2* tune. The SM cross sections and branching fractions are taken from Ref. |[20| . 

The samples of Drell-Yan events (qq — > Z/ 7 * — )■ , where £ is a lepton), and background 

samples used to represent the diphoton continuum and processes where one of the photon 
candidates arises from misidentified jet fragments, are the same as used in Ref. I1T2I . Simulated 
samples of Z — > e+e^, Z — > and Z — > events, used for comparison with data 

and to extract an energy scale and corrections for resolution of photon energies, are generated 
with MadGraph, SHERPA, and POWHEG 1401 , providing comparisons among the different 
generators. Simulated background samples are used for training of multivariate discriminants, 
and for defining selection and classification criteria. 

The diphoton continuum processes involving two prompt photons are simulated using SHERPA 
1.4.2 mil. The remaining processes where one of the photon candidates arises from misidenti¬ 
fied jet fragments are simulated using PYTHIA 6 alone. 

A comparison of unfolded data with results from models for ggH using the MG generators 
HRES Hill in, POWHEG and POWHEG+MINLO |42, and MADGRAPH5_aMC@NLO ||43 is pre¬ 
sented in Section 1^ 


4 Event selection and classification 

Trigger requirements, vertex determination, and kinematic criteria on photons are unchanged 
relative to those given in Ref. 112 1. The multivariate classifiers used to identify photons and to 
estimate mass resolution are also unchanged, but used in a different way. Instead of using a 
discriminant for photon identification as an input to the final diphoton kinematic discriminant, 
a requirement is set on the photon identification discriminant. Event classification, instead 
of being based on the output of the kinematic discriminant, is based on the estimated 
resolution. These differences are described in greater detail in the following section. 
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4 Event selection and classification 


4.1 Photon identification 

Photon candidates are required to be within the fiducial region of |?/| < 2.5, excluding the 
barrel-endcap transition region of 1.44 < \tj\ < 1.57, where photon reconstruction is not opti¬ 
mal. The transverse momenta of the two photons are required to satisfy the previously men¬ 
tioned conditions of /niry^ >1/3 and pj^/nijj >1/4. The use of pj thresholds scaled by 
niyry prevents the distortion of the low end of the nijj spectrum that results if a fixed thresh¬ 
old is used. Photons are also required to satisfy preselection criteria based on isolation and 
distributions in shower variables slightly more stringent than used in the trigger requirements. 

Three variables are calculated for each reconstructed candidate for a pp interaction vertex: the 
sum of the pj of the charged-particle tracks emerging from the vertex, and two variables that 
quantify the difference in the vector and scalar sums in pj between the diphoton system and 
the charged-particle tracks associated with the vertex. In addition, if either photon is associated 
with any charged-particle track identified as resulting from 7 —> e+e^ conversion, extrapola¬ 
tion of their trajectories is used to clarify the origin of the vertex of their production. These 
variables are used as inputs to a multivariate system based on a boosted decision tree (BDT) 
classifier to choose the reconstructed vertex to associate with the diphoton system. All BDTs 
are implemented using the TMVA Il46l framework. 

Another BDT is trained to separate prompt photons from photon candidates resulting from 
misidentification of jet fragments passing the preselection requirements. Inputs to the BDT 
are variables related to the lateral spread of the shower, and isolation energies reconstructed 
from scalar sums in pj of charged particles and Ej sums of photons in a cone with an opening 
angle AR = \/(Ap)'^ -|- {A(p)^ < 0.3 around the photon, computed using the particle-flow (PF) 
algorithm 1147111^ . 

The p and energy of the supercluster corresponding to the reconstructed photon are also in¬ 
cluded as input variables in the photon identification BDT. These variables are introduced to 
explicitly correlate the shower topology and isolation variables with p and pj. Furthermore, 
during the BDT training, the p and pj background distributions are reweighted to match the 
distributions in the signal. As a result, for a given requirement on the BDT output, efficiencies 
for photon identification are almost independent of p and pi, as can be seen in Fig. which 
provides less model-dependent efficiency corrections for comparison of data and MC expecta¬ 
tions at the particle level. A loose selection is applied on the BDT output for photons detected 
in the barrel region, while a tight selection is applied to endcap photons, with respective mean 
efficiencies of 95% and 90% for signal photons. 

Photon efficiencies in signal samples are corrected for the difference in efficiency between data 
and simulation as measured with Z —)■ e^e^ events, treating electrons as photons by reweight¬ 
ing the electron cluster variable Rg to that of the Rg distribution of signal photons. There is 
good agreement found in the BDT output between data and simulation in Z —)■ e+e^ and 
Z —events. 

4.2 Event classification using an estimator of mass resolution 

Measuring differential kinematic distributions implies that events cannot be classified accord¬ 
ing to the diphoton kinematic BDT used for the inclusive measurement of Ref. 11121 , as that 
would create a bias in the result. To improve the performance beyond the simple classification 
using the Rg variable in the reference sequential analysis, referred to as "cut-based analysis" in 
Ref. nil, an event categorization is introduced that is based on the estimated energy resolution. 

The photon energies are corrected using a multivariate regression technique for the contain- 
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Figure 1: Photon identification efficiencies for a Higgs boson with =125 GeV, as a function 
of photon pseudorapidity (left), and photon transverse momentum (right). 


ment of showers in clustered crystals, for shower losses of photons that convert in the material 
upstream of the calorimeter, and for effects from pileup, based on shower variables and vari¬ 
ables related to positions of photons in the detector studied in y+jets simulated events. The 
energy response to photons is parameterized through an extended form of the Crystal Ball 
function Il49ll with a Gaussian core and two power law contributions. The regression provides 
an estimate of the parameters of the function, and therefore a prediction for the distribution in 
the ratio of the uncorrected supercluster energy to the true energy The correction to the pho¬ 
ton energy is taken as the inverse of the most probable value of this distribution. The standard 
deviation of the Gaussian core provides an estimate of the uncertainty in the energy ((Je). 

We define the estimator of the diphoton mass resolution by: 



where £i and Ei are the corrected energies of the two photons from the regression, (Te^ and cte^ 
are the uncertainties in the photon energies, and Si and Sz are smearing terms depending on tj, 
Rg, and Ej, determined from Z —> e+e^ events, needed for the simulation to match the energy 
resolution in data. 

For the typical energy range of photons in this analysis, the relative energy resolution (rg / £ 
depends on the energy, which in turn introduces a dependence of the mass resolution on the 
value of ntyj. Therefore, a categorization based on the diphoton mass resolution introduces a 
distortion of the background mass spectrum. To obtain a smooth background description, we 
apply a transformation to the estimator (Te to decorrelate cte/E from the energy 

The value of c^e / £ depends on fj because of differences in material in front of the EGAL and the 
inherent properties of the EGAL. To ensure that the decorrelation is performed independent 
of the 7 / distribution of the training sample, in a first step a transformation is applied to make 
cte/E independent of £ and ?/. A 7 +jets MG sample is used to build the fully decorrelated 
variable, that covers a wide range in }j and pj. The decorrelation is performed by making a 
change of variable, replacing the probability distribution in cte / £ by its cumulative distribu¬ 
tion function cdf{crE/E). This function follows a uniform distribution f50 l and removes any 
correlation between (Te/E and {E,rj). Since the tj dependence is removed, this variable is not 
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suitable for estimating a per-event mass resolution, which is non-uniform in the detector. A 
second step is therefore introduced to restore the correlation of the energy resolution with pho¬ 
ton in a H —> 77 MC sample with niw = 123 GeV (statistically independent of the simulated 
events used to model the signal), and recover thereby a dependence of on rj for the photons 
of interest. In this step a new change of variable is performed, replacing the previous a^/E 
with the inverse of the cumulative distribution function oi a^/E{rj). The impact of the partic¬ 
ular Pj spectrum used for this step is only modest for the correlation of cr^ with t], which is 
in any case dominated by the material distribution in front of the ECAL and by calorimeter 
performance. The advantage of such a two-step procedure is that it offers a decorrelated cte/E 
variable that is uncorrelated with E and does not depend on the p distribution of the training 
sample. It provides the typical energy resolution at a given p of the detector. The final result 
can be interpreted as an estimator of the average energy resolution at a particular value of rj. 

The value in the estimator of the mass resolution, after decorrelation, is used to categorize 
events. The dm/niyy distribution in Z —> e+e^ data shown in Fig. [^indicates good agreement 
with the MC expectation over the whole range of crm/nijj. The double bump structure cor¬ 
responds mainly to events where (Fig. 2 a) both photons are in the central barrel {\f]\ < 1 . 0 , 
< 1%), or at least one of the photons is in the outer barrel (1.00 < \f]\ < 1.44, 
(Tm/mjy > 1%), and (Fig. 2b) one photon is in the barrel and one in the endcap, both hav¬ 
ing high values of the R 9 {cTm/mjj < 1.5%), and all the other events {am/myy > 1.5%). The 
class boundaries in CTmfinyy are optimized in MC samples simultaneously with the photon- 
identification working points desired to maximize signal significance. The photon identifi¬ 
cation efficiency is about 95% in each category of crm/niy^. The category with best energy 
resolution has (Tm/niyy < 0.79%, and is composed of events with both selected photons in 
the central barrel. Both the second (0.79 < am/nijy < 1.28%) and third categories (1.28 < 
(Tm/nijy < 1.83%) have at least one photon in the outer barrel or in the endcap. Events with 
CTm/niyy > 1.83% do not provide a noticeable improvement in sensitivity and are therefore re¬ 
jected. Classifying events in these three categories improves the analysis sensitivity by nearly 
10 %. 

It was verified in the simulation that the classification of events according to dm / after im¬ 
plementing the decorrelation procedure for (Te, does not produce a distortion of the background 
distribution in 

4.3 Jet identification 

Jets are reconstructed using particles identified by the PF algorithm, using the anti-kj IISTl algo¬ 
rithm with a distance parameter of 0.5. Jet energy corrections account in particular for pileup, 
and are obtained from simulation. They are calibrated with in situ measurements using the 
energy balance studied in dijet and 7 /Z+jet events Il52l . The jet momentum scale is found to be 
within 5-10% of the true jet momentum over the whole spectrum and detector acceptance. The 
jet energy resolution is typically 15% and 8 % at 10 and 100 GeV, respectively. Mean resolutions 
of 10% to 15% are observed in the respective regions of |?/| < 0.5 and 3 < \ri\ < 5. Jets are 
selected if they fail the pileup identification criteria l53l , and have pj > 25 GeV. The minimum 
distance between photons and jets is required to be AR( 7 ,j) = A^( 7 ,j )2 > 0.5, 

where A?^( 7 ,j) and A^( 7 ,j) are the pseudorapidity and azimuthal angle differences between 
photons and jets, to minimize photon energy depositions into jets. 
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Figure 2: Mass resolution estimator after the decorrelation procedure, in Z —> e+e^ 

events in data (dots) and simulated events (histogram) with their systematic uncertainties 
(shaded bands) for barrel-barrel events (left), all the other events (right). The ratio of data 
to MC predictions are shown below each panel, and the error bars on each point represent the 
statistical uncertainties of the data. 

4.4 Fiducial phase space and observables 

The generator-level fiducial volume is chosen to be close to that used in the selection of recon¬ 
structed events, and follows previous prescriptions: photons must have / niyy >1/3 and 
Pj^/niyy >1/4, with both photons within \t]\ < 2.5. The photons have to be isolated at the 
generator level, with Ejj <10 GeV, where i runs over all the other generator-level particles 
in a cone AR < 0.4 around the photons. This selection corresponds to a signal efficiency of 63% 
in ggH, and almost 60% considering other production mechanisms. We measure the kinematic 
observables using two-photon criteria, as well as requiring at least one or two jets. 

The transverse momentum pj'^ and absolute value of the rapidity of the Higgs boson are 
measured using the inclusive selection. Both and probe the production mechanism, 
while the photon helicity angle cos 0* in the Collins-Soper frame 1541 of the diphoton system 
and the difference in azimuth between the two photons are related to properties of the 
decaying particles. 

The number of jets N^ets, the transverse momentum pj of the jet with largest (leading) pj in the 
event, and the rapidity difference between the Higgs boson and the leading jet \y'^^ — yh| are 
defined after requiring at least one jet with pj > 25GeV to lie within \rj\ < 2.5. The difference 
in rapidities between the diphoton system and the leading jet provides a sensitive probe of any 
new contributions to the ggH process. 

Requiring at least two jets with pj > 25GeV and \r]\ < 4.7 provides the basis for defining the 
following observables: the dijet mass wjjj, the azimuthal angle difference of the two jets A(p^\ 
the difference in pseudorapidity Apij between the leading and subleading jet, the Zeppenfeld 
1551 variable — (//h -|- rj^'^')/2\, and the difference in azimuth between the Higgs boson and 
the dijet system A^'^'^')/ A requirement of |//^ | < 2.5 is applied in the single-jet selection, as this 
selection aims to probe primarily ggH process, while the two-jet selection has a requirement 
\i]^ < 4.7 since it is oriented toward VBF. 

The bin boundaries for each kinematic observable are optimized to achieve similar relative 
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5 Extraction of signal and unfolding of detector effects 


statistical uncertainties in the expected cross section in each bin, namely 60% for each bin in 
the inclusive observables, 70% in the one-jet observables, and more than 100% in the two-jet 
observables. The relative statistical uncertainty expected in the fiducial cross section is 30%. 


5 Extraction of signal and unfolding of detector effects 

The signal yield is extracted by fitting the m^y distribution using a signal model based on 
simulated events, and a background model determined in the fit to the data. The statistical 
methodology is similar to Ref. |[T2l , and for each observable the fit is performed simultaneously 
in all the bins. The reconstructed yields are corrected for detector effects by including the 
response matrix in the fit. 

5.1 Models for signal 

Signal models are constructed for each class of (Tm/niyy events, and for each production mech¬ 
anism, from a fit to the simulated niyy distribution, after applying the corrections determined 
from comparisons of data and simulation for Z —> e+e^ (also checked with Z —>• 
events, for wjh = 120,125, and 130 GeV. Mass distributions for the best and worst choices of 
diphoton vertex, corresponding to the highest and lowest scores in the vertexing BDT (121, are 
fitted separately. Good descriptions of the distributions can be achieved using sums of Gaus¬ 
sian functions, where the means are not required to be identical. As many as four contributing 
Gaussian functions are used, although in most cases two or three provide an acceptable fit. 
Models for intermediate values of mn are obtained by linear interpolation of the fitted param¬ 
eters. 

5.2 Statistical methodology 

After implementing the above-described selection requirements on photon candidates, a si¬ 
multaneous binned maximum likelihood fit is performed to the diphoton invariant mass dis¬ 
tributions in all the event classes over the range 100 < niyy < 180 GeV for each differential 
observable. The test statistic chosen to measure signal and background contributions in data is 
based on the profile likelihood ratio ||56j|57|. Systematic uncertainties are incorporated into the 
analysis via nuisance parameters and treated according to the frequentist paradigm. 

We use the same discrete profiling method to fit the background contribution (58l as used 
in extracting the main H —)■ 77 result (121 . The background is evaluated by fitting the niyy 
distribution in data, without reference to MG simulation. Thus the likelihood to be evaluated 
in a signal+background hypothesis is 

= C{data\si{iii,myy) -6 fi{myy)), ( 2 ) 

where jii is the signal strength (ratio of measured to expected yields) in the bin i of a differential 
distribution that is varied in the fit, Si{}ii,myy) represents the model for signal, and fi{myy) the 
fitted background functions. 

The choice of function used to fit the background in any particular event class is included 
as a discrete nuisance parameter in the formulation used to extract the result. Exponentials, 
power-law functions, polynomials in the Bernstein basis, and Laurent polynomials are used to 
represent f{myy). When fitting a signal+background hypothesis to the data, by minimizing 
the value of twice the negative logarithm of the likelihood, all functions in these families are 
tried, with a "penalty" term added to account for the number of free parameters in the fit. 
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The penalized likelihood function Cf for a single fixed background fitting function / is defined 
through 

-2\nCf = -2 \nCf + Nf, (3) 

where Cj is the "unpenalized" likelihood function, and Nj is the number of free parameters in 
/. The full set of denoted by fi, is determined by minimizing the likelihood ratio: 


L(^) = —2 In 


/:(data|^, Bfijfi) 

CidatalpJpJpY 


(4) 


where the numerator represents the maximum of £ as a function of ji, achieved for the best-fit 
values of the nuisance parameters 9^ = 9ji, and a particular background function = //;• The 
denominator corresponds to the global maximum of C, where jl = jl, 0^-; = 9p, and = fp. 
In each family, the number of degrees of freedom (number of exponentials, number of terms 
in the series, degree of the polynomial, etc.) is increased until no significant improvement 
(p value < 0.05 obtained from the F-distribution |^1) occurs in the likelihood between N+1 
and N degrees of freedom for the fit to data. 

For a given observable, the fit is performed simultaneously over all the bins, and the nuisance 
parameters are profiled in the fit. The signal mass is also considered a nuisance parameter 
and profiled for each observable. This choice is made to avoid using the same data twice, 
first to measure the signal mass, then to measure the kinematic distribution at this mass. As a 
consequence the differential cross section for each observable is evaluated at slightly different 
best fit values of the signal mass (see appendix A). As an example. Fig. shows the sum of the 
fit to the events of the three CTm/mjj classes in the fiducial phase space measurement, under the 
signal+background hypothesis (S+B), weighted by 5/(S+B) separately in each category 



(GeV) 

Figure 3: Sum of the signal+background (S+B) model fits to the events of the three (Tm/niyy 
classes in the fiducial phase space measurement, weighted by S/(S+B) separately in each cate¬ 
gory, together with the data binned as a function of niy-y. The 1 and 2 standard deviation bands 
of uncertainty (labeled as Ic and 20 ") shown for the background component include the uncer¬ 
tainty due to the choice of function and the uncertainty in the fitted parameters. The bottom 
panel shows the result after subtracting the background component. 
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5 Extraction of signal and unfolding of detector effects 


The uncertainty on the expected signal strength of the fiducial cross section is np = 0.32 for the 
present analysis, compared with the uncertainty on the expected signal strength obtained with 
the reference analysis of = 0.26 using 8 TeV data. The reference analysis Ifl2l classifies events 
using exclusive categories dedicated to measuring signal production in associated mechanisms, 
while the present analysis is performed inclusively 


5.3 Unfolding detector effects 

The measurement is performed simultaneously in all bins, together with the unfolding of the 
detector effects to the particle level. The same kind of procedure was used to extract the signal 
strength in untagged and dijet-tagged categories to measure the couplings of the Higgs bo¬ 
son to vector bosons and fermions 1121 . This procedure uses asymmetric uncertainties in the 
full likelihood instead of being limited to Gaussian uncertainties computed with a covariance 
matrix. 

The unfolding of reconstructed distributions is based on models of response matrices Kij in the 
three (Tm/mjj categories for each observable to be measured. The Kij are constructed in the 
simulation to give the probability of measuring a reconstructed event in bin j, given that it was 
generated in bin i. The models for contributions of signal are contained in each element of the 
response matrix, which is mostly diagonal, but can have non-negligible bin-by-bin migration 
resulting in off-diagonal contributions. 

The unfolding is performed using a maximum likelihood technique, adapted to combine the 
measurement in different categories. Regularization is not used, so as not to shift the best-fit 
value while artificially decreasing the uncertainties. This leads to minimizing the following 
conditional log-likelihood expression: 

H¥) = -2x:iog£(x:i^yH/NriNf“) (s) 

/ i 


where C is the log-likelihood expression in Eq. (2), is the unknown unfolded particle- 

level distribution, jii is the unknown signal strength at particle level, is the particle-level 
distribution in the simulated kinematic observable, and is the number of events in each 
bin of the measured distribution. The indices i refer to particle-level bins while j refer to 
reconstructed-level bins in the three categories. 

The expected reconstructed signal at detector level is given by the vector ]j = KijNf^^, where 
each entry of Kij corresponds therefore to a set of signal models, computed by interpolating 
the matrix between the generated mass points, weighted by the signal efficiency interpolated 
at the same mass. The matrix is a function of Wh, the model for the generated bin i and the 
reconstructed bin and category j, as well as all the nuisance parameters. Events falling out of 
the acceptance are taken into account by forming an extra bin. 


The maximum likelihood fit is performed simultaneously for the diphoton background and 
the signal strength in each generator bin, as described in Section 5.2 The out-of-acceptance 


bin is left fixed in the fit, because there are only N bins at the reconstructed level, where the 
detector is able to perform the measurement, and it is therefore not possible to determine an 
extra unknown at the particle level {N + 1 unknowns), outside of the detector acceptance. To 
restore the correct cross section normalization in the fiducial region, the generator distributions 
for the variables of interest (Nf^^ at the fitted mass point oth) are multiplied by the measured 
set of signal strengths {fij). 


The enhancement of the statistical uncertainties due to the presence of off-diagonal elements 
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in the response matrix is small for the inclusive observables, while it is non-negligible for the 
one-jet or two-jet observables. A negative number of events is measured in only one bin of the 
rapidity difference between the Higgs boson and the leading jet, in the last bin of the dijet mass 
distribution, and in two bins of the azimuthal difference between the Higgs boson and the dijet 
system. In each case, within the uncertainties, the result is compatible with zero. 

The model dependence introduced by the unfolding is checked using the following procedure. 
The same model used for the SM Higgs boson is kept in the unfolding matrix, leaving the ex¬ 
pected yield for ggH unchanged. The expected number of signal events arising from associated 
production mechanisms is altered by 50% in the fit. This change introduces a redistribution of 
the events in the CTm/myry categories relative to the nominal analysis. The change in the fiducial 
cross section is less than 5% of its statistical uncertainty. In general, the impact in each bin of 
the measured differential distributions can be up to 5% and 10% of the statistical uncertainty in 
the respective bins of the inclusive and jet observables. 


6 Systematic uncertainties 

Systematic uncertainties listed in this section are included in the likelihood as nuisance param¬ 
eters and are profiled during the minimization. Unless specified to the contrary, the sources of 
uncertainty refer to the individual quantity studied, and not to the final yield. The precision of 
the present measurement is however dominated by statistical uncertainties. 

The sources of uncertainty assigned to all events can be summarised as follows. The uncer¬ 
tainty in the integrated luminosity is estimated as described in Ref. 1601 , and amounts to 2.6% 
of the signal yield in the data. The uncertainty in the vertex-finding efficiency is taken from the 
difference observed between data and simulation in Z —> events, following removal of 

the muon tracks to mimic a diphoton event. A 1% uncertainty is added to account for the ac¬ 
tivity from charged-particle tracks in signal, estimated by changing the underlying event tunes 
in ggH events; another uncertainty of 0.2% accounts for the uncertainty in the distribution 
of signal events. The uncertainty in the trigger efficiency is extracted from Z —e+e~ events 
using a "tag-and-probe" technique 11611 . A rescaling in the Rg distribution is used to take into 
account the difference between electrons and photons, for a total uncertainty of 1 % assigned 
for this source. 

The following correspond to systematic uncertainties related to individual photons. The uncer¬ 
tainty in the energy scale of photons is assessed using simulated samples in which the amount 
of tracker material is increased uniformly by 10 % in the central barrel, where the material is 
known with best precision, and 20% out of this region. These values were chosen as upper 
limits on the additional material, as derived from the data. The resulting uncertainty in the 
photon energy ranges from 0.03% in the central ECAL barrel up to 0.3% in the outer endcap. 
Additional uncertainties of 0.015% are due to the modelling of the fraction of scintillation light 
reaching the photodetector, and from nonuniformities in the radiation-induced loss of trans¬ 
parency of the crystals. A small uncertainty of 0.05% is added to account for modelling of 
electromagnetic-showers in GEANT4 version 9.4.p03. 

Possible differences between MC simulation and data in the extrapolation of shower energies 
typical of electrons from Z ^ e+e^ decays to those typical of photons from H —> 77 decays, 
have been investigated with Z —)■ e+e^ and W —> ev data. The effect of differential nonlinearity 
in the measurement of photon energies has an effect of up to 0 . 1 % on the diphoton mass for 
niryr^ K, 125 GoV. 
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7 Comparison of data with theory 


The energy scale and resolution in data are measured with electrons from Z —> e+e^ decays. 
Systematic uncertainties in the method are estimated as a function of |) 7 | and Rg. The un¬ 
certainties range from 0.05% for unconverted photons in the central ECAL barrel, to 0.1% for 
converted photons in the outer endcaps of the ECAL. Einally, there is an overall uncertainty 
that accounts for possible mismodellmg of the Z —> e+e^ line shape in simulation. 

The uncertainties in the BDT discriminant for photon identification and in the estimate of pho¬ 
ton energy resolution are discussed together since they are studied in the same way. The 
dominant underlying cause of the observed differences between data and simulation is the 
simulation of the energy distribution in the shower. The combined contribution of the un¬ 
certainties in these two quantities dominates the experimental contribution to the system¬ 
atic uncertainty in signal strength. The agreement between data and simulation is examined 
when photon candidates are electron showers reconstructed as photons in Z —> e+e^ events, 
photons in Z —> events, and leading photons in preselected diphoton events where 

niryr^ > 160 GeV. A change of ±0.01 in the value of the photon identification discriminant, 
together with an uncertainty in the estimated photon energy resolution parameterized, respec¬ 
tively as a rescaling of the resolution estimate by ±5% and ±10% about its nominal value in the 
barrel and in the endcap, fully cover the differences observed in all three of these data samples. 

The uncertainty in the preselection efficiency is taken as the uncertainty in the data/MC scale 
factors measured using Z —> e+e^ events with a tag-and-probe technique. 

Jet observables are affected by systematic uncertainties arising from jet identification, jet energy 
scale, and resolution. Eor the jet observables, a systematic uncertainty of less than 1% in the 
impact of the algorithm used to reject jets from pileup is neglected. Eor jets within \r]\ <2.5, the 
energy scale uncertainty is «3% at 30 GeV, and decreases quickly as a function of the increasing 
jet Pt- The impact of this uncertainty in the cross section is 1-5%, increasing with the number 
of jets. Dijet observables for forward jets (up to |//1 < 4.7) have the worst energy resolution and 
a scale uncertainty of «4.5% at 30 GeV. The impact of these uncertainties on the cross section 
ranges from 5 to 11%, increasing in kinematical regions of observables where jet t] is large. 
Small contributions from corrections in jet energy resolution have an impact of less than 1% 
and are neglected. 

7 Comparison of data with theory 

7.1 Theoretical predictions 

The unfolded data are compared with the HRES Il42ll43l , POWHEG 127113011 , POWHEG+MINLO 1441 . 
and MADGRAPH5_aMC@NLO 145116^ MG generators for ggH production. 

The HRES parton-level generator corresponds to next-to-next-to-leading order (NNLO) accu¬ 
racy in perturbative QGD, with next-to-next-to-leading logarithm soft-gluon resummation; 
HRES v.2.3 assumes finite bottom and top quark masses, using respective values of ni}, = 
4.75 GeV and nit = 175 GeV. The renormalization scale ps and factorization scale pf are set 
to wh = 125 GeV, while the resummation scale is set to nin/l. The MSTW2008NNLO 1631 PDE 
is used for the central value, and its 68% confidence level eigenvectors for computing the uncer¬ 
tainty (following the LHAPDE l64l recipe). The dependence on scale is evaluated by changing 
independently both the renormalization and factorization scale up and down by a factor of two 
around the central value nin, and not considering simultaneous changes such as ps = niull 
and pf = 2?«h. Because HRES cannot be interfaced to a parton-shower program, no isolation is 
applied at the partonic level. A nonperturbative correction must be applied to the distributions 



7.1 Theoretical predictions 
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to correct for the efficiency loss due to isolation requirements in the presence of parton shower 
and underlying event. The nonperturbative correction is evaluated from the mean of the iso¬ 
lation efficiencies computed with powheg and MADGRAPH5_aMC@NLO (as described below), 
estimated to be 3.1% (up to 5% in some bins). The uncertainty is taken as half of the envelope, 
which is between 0.5% to 5%, depending on the kinematical region. 

The POWHEG parton-level generator implements NLO calculations 112711 interfaced to parton 
shower programs. Samples of events with a Higgs boson with nin = 125 GeV produced via 
ggH, assuming an infinite top quark mass, are generated and hadronized with PYTHIA 6.4. 
A sample of events with a Higgs boson produced in association with just a single jet, called 
POWHEG HJ, is also generated with POWHEG+MINLO HMJ- This sample has NLO accuracy for 
0-jet and 1-jet production, while it is only leading-order (LO) for 2-jet final states. Both samples 
set the damping factor hfact in POWHEG at 100 GeV to reproduce the predicted pj distribution 
of the Higgs boson from HqT Il35ll36l . This factor minimizes emission of extra jets beyond those 
in the matrix element in the limit of large pj, and enhances contribution from the POWHEG 
Sudakov form factor as px approaches 0. The GTIO PDF and PYTHIA 6 tune 22* are used in the 
calculation. Theoretical uncertainties are computed in the same way as described for hres. 

The MADGRAPH5_aMC@NLO matrix element generator is capable of generating LO and NLO 
processes 1162| . The ggH process is generated using the NLO Higgs characterization model 
ll65l , with effective coupling of the Higgs boson to gluons in the infinite top quark mass limit. 
Gluon fusion is generated with 0,1, or 2 additional jets at NLO in the Born matrix element, and 
combined using FxFx merging Il 66 l . Samples are generated using the GTIO PDF, and show¬ 
ered using PYTHIA 8.185 167 1 with the 4G tune 168 1. A nominal merging scale of 30 GeV is 
used for the additional jet multiplicities. The effect of changing the merging scale from 20 to 
60 GeV is very small compared to uncertainties in scale and in the choice of PDF. Uncertainties 
from renormalization and factorization scales are evaluated in the same way as with HRES and 
POWHEG. 

Theoretical predictions for associated production mechanisms are computed with the follow¬ 
ing generators. POWHEG interfaced with PYTHIA 6 is used for VBF, while standalone PYTHIA 
6 is used for VH and tfH. In the following, the notation XH refers to the sum of VBF, VH 
and tfH predictions for these generators. Each of the ggH predictions for HRES, POWHEG, 
MADGRAPH5_aMC@NLO, POWHEG+MINLO, and XH processes are normalized to the total cross 
sections from Ref. Il20l . 

Along with the SM predictions, the following alternative models are considered. Spin 2+ mini¬ 
mal model (graviton-like) initiated through two production mechanisms: ggH and qq annihil- 
iation, based on the JHUGEN generator [2111691 and normalized to the total SM cross section. 
The main changes relative to the SM are expected in the inclusive Gollins-Soper cos 6* angu¬ 
lar distribution in the 77 rest frame, which provides maximum information on the spin of 
the 77 system. The two spin-2 samples are compared to the data in the cos0* observable. 
Anomalous couplings parametrized with the Ow dimension -6 operator in linearly realized ef¬ 
fective field theory [70| are also considered, and implemented through the Universal FeynRules 
Output llTIII in MadGraph 5. The operator is related to the anomalous triple gauge cou¬ 
pling parameter A^f [ 721 . The values of the Wilson coefficients are Fw = —5 x lO^^GeV^"^ 
and Fw = +5 X 10^^ GeV^"^, both corresponding to Agf = 0.21, a value approximately five 
times the size of the limits set by LHG diboson measurements Il73] - [76l . Both values modify 
the kinematic distributions of the Higgs boson in the VBF process toward larger pj, and the 
Fyj = —5 X 10^^ GeV^^ value also increases the VBF cross section by approximately a factor of 
3. 
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7 Comparison of data with theory 


All the predictions are generated at mu = 125 GeV. For each observable, a correction fac¬ 
tor is applied in each bin of the differential cross section to correct for the mass difference of 
the generated sample relative to the measured in data. The correction is computed with 
POWHEG+PYTHIA 8 for both the ggH and VBF processes, and with PYTHIA 6 for VH and tfH 
processes. It amounts to less than 1% for all bins, and integrates to a 0.8% effect in the fiducial 
cross section. 

7.2 Results 

The fiducial cross section, inclusive in the number of jets, is measured to be: 

%bs = 32+|o (stat )+3 (syst) fb, 

where the uncertainties reflect statistical and systematic contributions added in quadrature. 
This can be compared with the following SM predictions: 

Ghres+xh = Sll^fb, 

GpoWHEG+XH = ^'^-5 ft), 

PMADGRAPH5_flMC@NLO+XH = 30^5 ft- 

Uncertainties in the predicted cross sections include contributions from renormalization and 
factorization scales, choice of PDF and branching fraction. The hres also includes uncertainties 
from nonperturbative corrections. 

The observed fiducial cross section agrees with the predicted values. The measurement pre¬ 
cision is dominated by statistical uncertainties. The relative systematic uncertainty of 9% is 
almost negligible relative to the statistical uncertainty of 30%. The experimental uncertainty is 
larger than the theoretical one by a factor up to about two. The ratio of the measured cross sec¬ 
tion to the predictions for POWHEG+XH is in good agreement with the signal strength observed 
in Ref. 1121. 

The measured differential cross sections observed in data, given for each bin by are 

compared with predictions for inclusive production in Fig.|^ and for jet observables in Figs. 
and 1^ The total theoretical uncertainty included in these comparisons is computed by adding 
in an uncorrelated way the uncertainties in the choice of PDF, renormalization and factor¬ 
ization scale, and the branching fraction. The uncertainties in the ggH mechanism, PDF 
choice, and the renormalization and factorization scales are computed with HRES, POWHEG, 
POWHEG+MINLO and MADGRAPH5_aMC@NLO, as described above, while the branching frac¬ 
tion for all production mechanisms, as well as the scale and PDF for the associated produc¬ 
tion mechanisms are taken from Ref. [20]. Distributions for inclusive observables computed 
with HRES, POWHEG, and MADGRAPH5_aMC@NLO, and jet-related observables computed with 
POWHEG, POWHEG+MINLO and MADGRAPH5_aMC@NLO, including latest higher-order correc¬ 
tions, are compatible within their uncertainties, and also compatible with the data. 

Figure]^ (upper left) shows the pj distribution of the Higgs boson, which is sensitive to higher- 
order corrections in perturbative QCD. Figure (upper right) shows the absolute rapidity dis¬ 
tribution of the Higgs boson, which is sensitive to the proton PDF, as well as to the production 
mechanism. Figure (lower left) shows the distribution. Figure (lower right) displays 
the cos 9* distribution, which is sensitive to the spin of the Higgs boson. The two spin-2 sam¬ 
ples indicate deviations relative to the SM predictions. As in the case of Ref. fl2l . the data do 
not have sufficient sensitivity to discriminate between spin-2 and spin-0 hypotheses. 
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Figure 1^ (upper left) shows the pj distribution for the leading jet, which is sensitive to higher- 
order QCD effects, and Fig. (upper right) shows the rapidity difference between the Higgs bo¬ 
son and the leading jet. The distribution in the number of jets is displayed in Fig.|^(lower left). 
The last bin gives the cross section for signal events with at least three jets. The distribution of 
the dijet mass shown in Fig. (lower right) is sensitive to the VBF production mechanism, and 
especially to a possible anomalous electroweak production in the tail of the distribution. Large 
contributions from new processes modifying triple gauge couplings would be detected as ex¬ 
cesses either in p^^, pj or in mjj distributions. The distributions in data are compatible with 
expectations from the SM within their uncertainties. Figure]^ (upper left) shows the difference 
in azimuthal angle between the two jets of highest pj. Figures (upper right), (lower 
left), and 1 ^ (lower right) show, respectively, the distribution of the rapidity difference between 
the two jets, the Zeppenfeld variable, and the azimuthal difference between the Higgs boson 
and the dijet system. These angular variables are sensitive to the VBF topology, and large con¬ 
tributions from anomalous couplings are not observed in data. The distributions in data are 
compatible with SM predictions within their statistical, systematic, and theoretical uncertain¬ 
ties. 


8 Summary 

A measurement was carried out of differential cross sections as a function of kinematic observ¬ 
ables in the H ^ 77 decay channel, using data collected by the CMS experiment at = 8 TeV, 
corresponding to an integrated luminosity of 19.7fb~^. The measurement was performed for 
events with two isolated photons in the kinematic range pj^ >1/3 and pj^/rtijj >1/4, 
with photon pseudorapidities within |j/1 < 2.5. Photon identification was chosen to reduce the 
dependence of the measurement on the kinematics of the signal. Event classification relied on 
an estimator of diphoton mass resolution. The signal extraction and the unfolding of experi¬ 
mental resolution were performed simultaneously in all bins of the chosen observables. In this 
kinematic range, the fiducial cross section was measured to be 32 ± 10 fb. 

The differential cross section of the Higgs boson was measured, inclusively in the number of 
jets, as a function of its transverse momentum p'Jp, its rapidity the Collins-Soper angu¬ 
lar variable cos 6*, the difference in azimuthal angle between the two photons and the 

number of associated jets Njets- The transverse momentum of the leading jet pj, and the dif¬ 
ference in rapidity between the Higgs boson and the leading jet | were determined in 

events with at least one accompanying jet. In events with at least two jets, measurements were 
made of the dijet mass mjj, the azimuth between the two jets the pseudorapidity difference 
between the two jets A/yii, the Zeppenfeld variable \ + rf>^)/2\, and the azimuthal an¬ 

gle between the Higgs boson and the dijet system A^'^'^'U. The differential cross sections were 
compared with several SM and beyond SM calculations, and found to be compatible with the 
SM predictions within statistical, systematic, and theoretical uncertainties. With more data, 
anomalous couplings of the Higgs boson could be measured from its differential distributions. 
It would be possible to discriminate among different SM predictions for Higgs boson produc¬ 
tion those that are providing the best description. 


Acknowledgments 

We congratulate our colleagues in the CERN accelerator departments for the excellent perfor¬ 
mance of the LHC and thank the technical and administrative staffs at CERN and at other 



16 


8 Summary 


> 

O 

O 

]q 


Q. 

■D 

'B 

■D 


10 " 


10 


10 '^ 


10 ' 


19.7fb'^ (8TeV) 


CMS 


J' 


* Data 

as^^(H^ 7 Y) from LHC Higgs XS WG 
H aMC@NLO + XH 
'^1$ HRES + XH 
POWHEG + XH 

- POWHEG + XH, FW*5E"® Gev’ 

- POWHEG + XH, FW«-5E"® GeV"' 

I-1 XH = VBF+VH+tfH 








“L 




50 


100 


150 




-IS 




200 




"D 

ro 


: 140 


":&120- 
<1 

|iooh 


X 

X 


19.7fb~' (8TeV) 


CMS 


3r 

2t- 

1 

0 


* Data 

a 3 ^(H^YY) from LHC Higgs XS WG 
H aMC@NLO + XH 
HRES + XH 
POWHEG + XH 

- POWHEG + XH, FW*5E"® GeV"^ 

- POWHEG + XH, FW.-5E"® GeV 

I-1 XH = VBF+VH+tfH 


r«i 1 


0 


m 






2 3 

Ajy I 


19.7tb ' (STeV) 


A(t)" [rad] 



Figure 4: The H —> 77 differential cross section for inclusive events as a function of (upper 
left) pj'^, (upper right) (lower left) Acp'^'^, and (lower right) |cos0*|. All the SM contri¬ 

butions are normalized to their cross section from Ref. [20|. Theoretical uncertainties in the 
renormalization and factorization scales, PDF, and branching fraction are added in quadrature. 
The error bars on data points reflect both statistical and systematic uncertainties. The last bin 
of distribution sums the events above 200 GeV. For each graph, the bottom panel shows 
the ratio of data to theoretical predictions from the POWHEG generator. 
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Figure 5: The H —)■ 77 differential cross section for H+jets events as a function of (upper left) 
Njets/ (upper right) Pj, (lower left) — y^^\, with jets within \rj\ < 2.5, and (lower right) 
njjj with jets within \rj\ < 4.7. All the SM contributions are normalized to their cross section 
from Ref. [1 20 ^^1. Theoretical uncertainties in the renormalization and factorization scales, PDF, 
and branching fraction are added in quadrature. The error bars on data points reflect both 
statistical and systematic uncertainties. In each distribution, the last bin corresponds to the 
sum over the events beyond the bins shown in the figure. For each graph, the bottom panel 
shows the ratio of data to theoretical predictions from the POWHEG generator. 
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Figure 6 : The H —> 77 differential cross section for H+2 jets events, with jets within \ t]\ < 4.7, as 
a function of (upper left) (upper right) (lower left) Zeppenfeld variable, and (lower 
right) A^'^'T'U. All the SM contributions are normalized to their cross section from Ref. Il20l . 
Theoretical uncertainties in the renormalization and factorization scales, PDF, and branching 
fraction are added in quadrature. The error bars on data points reflect both statistical and sys¬ 
tematic uncertainties. In the Zeppenfeld distribution, the last bin sums the events with values 
above 1.2. For each graph, the bottom panel shows the ratio of data to theoretical predictions 
from the POWHEG generator. 





































































































19 


CMS institutes for their contributions to the success of the CMS effort. In addition, we grate¬ 
fully acknowledge the computing centers and personnel of the Worldwide LHC Computing 
Grid for delivering so effectively the computing infrastructure essential to our analyses. Fi¬ 
nally, we acknowledge the enduring support for the construction and operation of the LHC 
and the CMS detector provided by the following funding agencies: the Austrian Federal Min¬ 
istry of Science, Research and Economy and the Austrian Science Fund; the Belgian Fonds de 
la Recherche Scientifique, and Fonds voor Wetenschappelijk Onderzoek; the Brazilian Fund¬ 
ing Agencies (CNPq, CAPES, EAPERJ, and EAPESP); the Bulgarian Ministry of Education and 
Science; CERN; the Chinese Academy of Sciences, Ministry of Science and Technology, and Na¬ 
tional Natural Science Eoundation of China; the Colombian Eunding Agency (COLCIENCIAS); 
the Croatian Ministry of Science, Education and Sport, and the Croatian Science Eoundation; 
the Research Promotion Foundation, Cyprus; the Ministry of Education and Research, Esto¬ 
nian Research Council via IUT23-4 and IUT23-6 and European Regional Development Eund, 
Estonia; the Academy of Einland, Einnish Ministry of Education and Culture, and Helsinki 
Institute of Physics; the Institut National de Physique Nucleaire et de Physique des Partic- 
ules / CNRS, and Commissariat a TEnergie Atomique et aux Energies Alternatives / CEA, 
Erance; the Bundesministerium fiir Bildung und Eorschung, Deutsche Eorschungsgemeinschaft, 
and Helmholtz-Gemeinschaft Deutscher Eorschungszentren, Germany; the General Secretariat 
for Research and Technology, Greece; the National Scientific Research Eoundation, and Na¬ 
tional Innovation Office, Hungary; the Department of Atomic Energy and the Department 
of Science and Technology, India; the Institute for Studies in Theoretical Physics and Mathe¬ 
matics, Iran; the Science Foundation, Ireland; the Istituto Nazionale di Fisica Nucleare, Italy; 
the Ministry of Science, IGT and Future Plarming, and National Research Foundation (NRF), 
Republic of Korea; the Lithuanian Academy of Sciences; the Ministry of Education, and Uni¬ 
versity of Malaya (Malaysia); the Mexican Eunding Agencies (GINVESTAV, GONAGYT, SEP, 
and UASLP-FAI); the Ministry of Business, Innovation and Employment, New Zealand; the 
Pakistan Atomic Energy Gommission; the Ministry of Science and Higher Education and the 
National Science Gentre, Poland; the Fundagao para a Giencia e a Tecnologia, Portugal; JINR, 
Dubna; the Ministry of Education and Science of the Russian Eederation, the Eederal Agency of 
Atomic Energy of the Russian Eederation, Russian Academy of Sciences, and the Russian Eoun¬ 
dation for Basic Research; the Ministry of Education, Science and Technological Development 
of Serbia; the Secretaria de Estado de Investigacion, Desarrollo e Irmovacion and Programa 
Gonsolider-Ingenio 2010, Spain; the Swiss Funding Agencies (ETH Board, ETH Zurich, PSI, 
SNF, UniZH, Ganton Zurich, and SER); the Ministry of Science and Technology, Taipei; the 
Thailand Genter of Excellence in Physics, the Institute for the Promotion of Teaching Science 
and Technology of Thailand, Special Task Force for Activating Research and the National Sci¬ 
ence and Technology Development Agency of Thailand; the Scientific and Technical Research 
Gouncil of Turkey, and Turkish Atomic Energy Authority; the National Academy of Sciences 
of Ukraine, and State Eund for Eundamental Researches, Ukraine; the Science and Technology 
Eacilities Gouncil, UK; the US Department of Energy, and the US National Science Eoundation. 

Individuals have received support from the Marie-Gurie program and the European Research 
Gouncil and EPLANET (European Union); the Leventis Eoundation; the A. P. Sloan Eounda¬ 
tion; the Alexander von Humboldt Eoundation; the Belgian Eederal Science Policy Office; the 
Eonds pour la Eormation a la Recherche dans ITndustrie et dans TAgriculture (ERIA-Belgium); 
the Agentschap voor Irmovatie door Wetenschap en Technologie (IWT-Belgium); the Ministry 
of Education, Youth and Sports (MEYS) of the Gzech Republic; the Gouncil of Science and In¬ 
dustrial Research, India; the HOMING PLUS program of the Eoundation for Polish Science, 
cofmanced from European Union, Regional Development Eund; the OPUS program of the Na¬ 
tional Science Genter (Poland); the Gompagnia di San Paolo (Torino); the Gonsorzio per la 



20 


References 


Fisica (Trieste); MIUR project 20108T4XTM (Italy); the Thalis and Aristeia programs cofinanced 
by EU-ESF and the Greek NSRF; the National Priorities Research Program by Qatar National 
Research Fund; the Rachadapisek Sompot Fund for Postdoctoral Fellowship, Chulalongkorn 
University (Thailand); and the Welch Foundation, contract C-1845. 


References 

[1] ATLAS Collaboration, "Observation of a new particle in the search for the Standard 
Model Higgs boson with the ATLAS detector at the LHC", Phys. Lett. B 716 (2012) 1, 

doi:10.1016/j.physletb.2 012.08.020, arXiv:12 07.7214 

[2] CMS Collaboration, "Observation of a new boson at a mass of 125 GeV with the CMS 
experiment at the LHC", Phys. Lett. B 716 (2012) 30, 

doi:10.1016/j.physletb.2012.08.021 

[3] CMS Collaboration, "Observation of a new boson with mass near 125 GeV in pp 
collisions at ^/s = 7 and 8 TeV", JHEP 06 (2013) 081, 

doi:10.1007/JHEP0 6(2 013) 081, arXiv:1303.4571 

[4] S. L. Glashow, "Partial-symmetries of weak interactions", Nucl. Phys. 22(1961)579, 

doi:10.1016/0029-5582(61)90469-2 

[5] S. Weinberg, "A Model of Leptons", Phys. Rev. Lett. 19 (1967) 1264, 

doi:10.1103/PhysRevLett.19.1264 

[6] A. Salam, "Weak and electromagnetic interactions", in Elementary particle physics: 
relativistic groups and analyticity, N. Svartholm, ed., p. 367. Almqvist & Wiskell, 1968. 
Proceedings of the eighth Nobel symposium. 

[7] E Englert and R. Brout, "Broken Symmetry and the Mass of Gauge Vector Mesons", 
Phys. Rev. Lett. 13 (1964) 321, doi : 10.1103/PhysRevLett .13.321 

[8] P. W. Higgs, "Broken Symmetries and the Masses of Gauge Bosons", Phys. Rev. Lett. 13 
(1964) 508, doi : 10.1103/PhysRevLett. 13.508, 

[9] G. S. Guralnik, C. R. Hagen, and T. W. B. Kibble, "Global Conservation Laws and 
Massless Particles", Phys. Rev. Lett. 13 (1964) 585, 

doi:10.1103/PhysRevLett.13.585 

[10] P. W. Higgs, "Spontaneous Symmetry Breakdown without Massless Bosons", Phys. Rev. 
145 (1966) 1156, doi : 10.1103/PhysRev. 145.1156 

[11] T. W. B. Kibble, "Symmetry Breaking in Non-Abelian Gauge Theories", Phys. Rev. 155 
(1967) 1554, doi : 10.1103/PhysRev. 155.1554 

[12] CMS Collaboration, "Observation of the diphoton decay of the Higgs boson and 
measurement of its properties", Eur. Phys. }. C 74 (2014) 3076, 

doi:10.1140/epjc/sl0052-014-3076-z,arXiv:1407.0558 

[13] ATLAS Collaboration, "Measurement of Higgs boson production in the diphoton decay 
channel in pp collisions at center-of-mass energies of 7 and 8 TeV with the ATLAS 
detector", Phys. Rev. D 90 (2014) 112015, doi : 10.1103/PhysRevD .90.112015, 
arXiv:1408.7084 




References 


21 


[14] CMS Collaboration, "Precise determination of the mass of the Higgs boson and tests of 
compatibility of its couplings with the standard model predictions using proton 
collisions at 7 and 8 TeV", Eur. Phys. J. C 75 (2015) 212, 

doi:10.1140/epjc/sl0052-015-3351-7, arXiv:1412.8 662 

[15] ATLAS Collaboration, "Evidence for the spin-0 nature of the Higgs boson using ATLAS 
data", Phys. Lett. B 726 (2013) 120, doi : 10.1016/j .physletb. 2013.08.026, 
arXiv:1307.1432 

[16] CMS Collaboration, "Constraints on the spin-parity and anomalous HVV couplings of 
the Higgs boson in proton collisions at 7 and 8 TeV", Phys. Rev. D 92 (2015), no. 1, 
012004, doi:10.1103/PhysRevD.92.012004, arXiv:1411.3441 

[17] ATLAS Collaboration, "Measurements of fiducial and differential cross sections for 
Higgs boson production in the diphoton decay channel at ^/s = 8 TeV with ATLAS", 
]HEP 09 (2014) 112, doi : 10.1007/JHEP09 (2014) 112, arXiv: 1407.4222 

[18] ATLAS Collaboration, "Fiducial and differential cross sections of Higgs boson 
production measured in the four-lepton decay charmel in pp collisions at ^/s = 8 TeV 
with the ATLAS detector", Phys. Lett. B 738 (2014) 234, 

doi:10.1016/j.physletb.2 014.09.054, arXiv:1408.3226 

[19] ATLAS Collaboration, "Measurements of the Total and Differential Higgs Boson 
Production Cross Sections Combining the H —> 77 and H —> ZZ* —)■ 4£ Decay Charmels 
at -\/s = 8 TeV with the ATLAS Detector", (2015). arXiv :1504.05833, Submitted to 
PRL. 

[20] LHC Higgs Cross Section Working Group, "Handbook of LHC Higgs cross sections: 3. 
Higgs Properties", CERN Report CERN-2013-004,2013. 

doi:10.5170/CERN-2013-004,arXiv:1307.1347, 

[21] Y. Gao et al., "Spin determination of single-produced resonances at hadron colliders", 
Phys. Rev. D 81 (2010) 075022, doi : 10.1103/PhysRevD .81.075022, 

arXiv:1001.3396 

[22] G. Grojean, E. Salvioni, M. Schlaffer, and A. Weiler, "Very boosted Higgs in gluon fusion", 
]HEP 05 (2014) 022, doi : 10.1007/JHEP05 (2014) 022, arXiv: 1312.3317 

[23] GMS Gollaboration, "The GMS experiment at the GERN LHG", JINST 3 (2008) S08004, 

doi:10.1088/1748-0221/3/08/S08004 

[24] GMS Gollaboration, "Energy calibration and resolution of the GMS electromagnetic 
calorimeter in pp collisions at = 7 TeV", JINST 8 (2013) P09009, 

doi:10.1088/1748-0221/8/09/P09009,arXiv:1306.2016 

[25] GMS Gollaboration, "Performance of photon reconstruction and identification with the 
GMS detector in proton-proton collisions at -\/s = 8 TeV", JINST 10 (2015) P08010, 

doi:10.1088/1748-0221/10/08/P08010, arXiv:1502.02 7 02 

[26] GEANT4 Gollaboration, "GEANT4—a simulation toolkit", Nucl. Instrum. Meth. A 506 
(2003) 250, doi : 10.1016/S0168-90 02 (03) 01368-8, 

[27] P Nason, "A new method for combining NLO QGD with shower Monte Garlo 
algorithms", JHEP 11 (2004) 040, doi : 10.1088/1126-6708/2004/11/040, 
arXiv:hep-ph/0409146 





22 


References 


[28] S. Frixione, P. Nason, and C. Oleari, "Matching NLO QCD computations with parton 
shower simulations: the POWHEG method", JHEP 11 (2007) 070, 

doi:10.1088/1126-6708/2007/11/07 0, arXiv:070 9.20 92 

[29] S. Alioli, P Nason, C. Oleari, and E. Re, "A general framework for implementing NLO 
calculations in shower Monte Carlo programs: the POWHEG BOX", JHEP 06 (2010) 043, 

doi:10.1007/JHEP0 6(2 010) 043, arXiv:1002.2 581 

[30] S. Alioli, P Nason, C. Oleari, and E. Re, "NLO Higgs boson production via gluon fusion 
matched with shower in POWHEG", JHEP 04 (2009) 002, 

doi:10.1088/1126-6708/2009/04/002, arXiv:0812.057 8 

[31] P Nason and C. Oleari, "NLO Higgs boson production via vector-boson fusion matched 
with shower in POWHEG", JHEP 02 (2010) 037, doi : 10.1007/JHEP02 (2010) 037, 
arXiv:0911.5299 

[32] T. Sjdstrand, S. Mrenna, and P Z. Skands, "PYTHIA 6.4 physics and manual", JHEP 05 
(2006) 026, doi : 10.1088/112 6-67 08/2 006/05/02 6, arXiv :hep-ph/0 603175, 

[33] H.-L. Lai et al., "New parton distributions for collider physics", Phys. Rev. D 82 (2010) 
074024, doi : 10.1103/PhysRevD . 82.074 024, arXiv: 1007.22 41 

[34] LHC Higgs Cross Section Working Group, "Handbook of LHC Higgs cross sections: 2. 
Differential Distributions", CERN Report CERN-2012-002, 2012. 

doi:10.5170/CERN-2012-002,arXiv:1201.3084, 

[35] G. Bozzi, S. Catani, D. de Elorian, and M. Grazzmi, "Transverse-momentum 
resummation and the spectrum of the Higgs boson at the EHC", Nucl. Phys. B 737 (2006) 
73, doi:10.1016/j.nuclphysb.2005.12.022,arXiv:hep-ph/05080 68 

[36] D. de Elorian, G. Eerrera, M. Grazzini, and D. Tommasini, "Transverse-momentum 
resummation: Higgs boson production at the Tevatron and the LHC", JHEP 11 (2011) 
064, doi:10.1007/JHEP11(2011)064,arXiv:1109.2109 

[37] R. Eield, "Early LHC Underlying Event Data - Eindings and Surprises", in 22nd Hadron 
Collider Physics Symposium (HCP 2010), W. Trischuk, ed. Toronto, 2010. 

arXiv:1010.3558 

[38] J. Pumplin et al., "New generation of parton distributions with uncertainties from global 
QCD analysis", JHEP 07 (2002) 012, doi : 10.1088/1126-6708/2002/07/012, 
arXiv:hep-ph/0201195 

[39] L. J. Dixon and M. S. Siu, "Resonance-Continuum Interference in the Diphoton Higgs 
Signal at the EHC", Phys. Rev. Lett. 90 (2003) 252001, 

doi:10.1103/PhysRevLett.90.252001, arXiv:hep-ph/0302233 

[40] E. Re, "NLO corrections merged with parton showers for Z+2 jets production using the 
POWHEG method", JHEP 10 (2012) 031, doi : 10.1007/JHEPIO (2 012) 031, 
arXiv:1204.5433 

[41] T. Gleisberg et al., "Event generation with SHERPA 1.1", JHEP 02 (2009) 007, 

doi:10.1088/1126-6708/2009/02/007, arXiv:0811.4 622 





References 


23 


[42] D. de Florian, G. Ferrera, M. Grazzini, and D. Tommasini, "Higgs boson production at 
the LHG: transverse momentum resummation effects in the H —> 77 , H —> WW —> Ivlv 
and H —> ZZ —)■ 4Z decay modes", JHEP 06 (2012) 132, 

doi:10.1007/JHEP0 6(2 012)132, arXiv:12 03.6321 

[43] M. Grazzini and H. Sargsyan, "Heavy-quark mass effects in Higgs boson production at 
the LHG", JHEP 09 (2013) 129, doi : 10.1007/JHEP09 (2013) 129, 

arXiv:1306.4581 

[44] K. Hamilton, P. Nason, and G. Zanderighi, "MINLO: multi-scale improved NLO", JHEP 
10 (2012) 155, doi : 10.1007/JHEP10 (2012)155, arXiv : 1206.3572 

[45] S. Frixione and B. R. Webber, "Matching NLO QGD computations and parton shower 
simulations", JHEP 06 (2002) 029, doi : 10.1088/1126-6708/2002/06/029, 
arXiv:hep-ph/0204244 

[46] H. Voss, A. Hocker, J. Stelzer, and F. Tegenfeldt, "TMVA: Toolkit for Multivariate Data 
Analysis with ROOT", in Xlth International Workshop on Advanced Computing and Analysis 
Techniques in Physics Research (ACAT), p. 40. 2007. arXiv : physics/0703039, 

[47] GMS Gollaboration, "Particle-Flow Event Reconstruction in GMS and Performance for 
Jets, Taus, and MET", GMS Physics Analysis Summary GMS-PAS-PET-09-001, 2009. 

[48] GMS Gollaboration, "Gommissioning of the Particle-flow Event Reconstruction with the 
first LHG collisions recorded in the GMS detector", GMS Physics Analysis Summary 
GMS-PAS-PET-10-001, 2010. 

[49] M. J. Oreglia, "A study of the reactions tp' —>• JJip" PhD thesis, Stanford University, 

1980. SLAG Report SLAG-R-236, see Appendix D. 

[50] G. Gowan, "Statistical Data Analysis". Oxford science publications. Glarendon Press, 
1998. 

[51] M. Gacciari, G. P. Salam, and G. Soyez, "The anti-kf jet clustering algorithm", JHEP 04 
(2008) 063, doi : 10.1088/112 6-6708/2008/04/0 63, arXiv: 0802.118 9 

[52] GMS Gollaboration, "Determination of jet energy calibration and transverse momentum 
resolution in GMS", JINST 6 (2011) P11002, 

doi:10.1088/1748-0221/6/11/P11002,arXiv:1107.4277 

[53] GMS Gollaboration, "Pileup Jet Identification", GMS Physics Analysis Summary 
GMS-PAS-JME-13-005,2013. 

[54] J. G. Gollins and D. E. Soper, "Angular distribution of dileptons in high-energy hadron 
collisions", Phys. Rev. D 16 (1977) 2219, doi : 10.1103/PhysRevD .16.2219 

[55] D. E. Rainwater, R. Szalapski, and D. Zeppenfeld, "Probing color-singlet exchange in 
Z -|- 2-jet events at the GERN LHG", Phys. Rev. D 54 (1996) 6680, 

doi:10.1103/PhysRevD.54.6 68 0, arXiv:hep-ph/9 6054 4 4 

[56] ATLAS and GMS Gollaborations, LHG Higgs Gombmation Group, "Procedure for the 
LHG Higgs boson search combination in Summer 2011", Technical Report 
ATL-PHYS-PUB 2011-11, GMS NOTE 2011/005, GERN, 2011. 





24 


References 


[57] CMS Collaboration, "Combined results of searches for the standard model Higgs boson 
in pp collisions at y/s = 1 TeV", Phys. Lett. B 710 (2012) 26, 

doi:10.1016/j.physletb.2 012.02.0 64, arXiv:12 02.14 88 

[58] P. D. Dauncey, M. Kenzie, N. Wardle, and G. J. Davies, "Handling uncertainties in 
background shapes: the discrete profiling method", JINST 10 (2015) P04015, 

doi:10.1088/1748-0221/10/04/P04015,arXiv:1408.6865 

[59] E. W. Weisstein, "F-Distribution". 

http : //mathworld. wolfram. com/F-Distribution . html, 2014. From 
Math World - A Wolfram Web Resource. 


[60] CMS Collaboration, "CMS Fuminosity Based on Pixel Cluster Counting - Summer 2013 
Update", CMS Physics Analysis Summary CMS-PAS-FUM-13-001, 2013. 

[61] CMS Collaboration, "Measurement of the Inclusive W and Z Production Cross Sections 
in pp Collisions at ^/s = 7 TeV", JHEP 10 (2011) 132, 

doi:10.1007/JHEP10(2 011)132, arXiv:1107.478 9 

[62] J. Alwall et al., "The automated computation of tree-level and next-to-leading order 
differential cross sections, and their matching to parton shower simulations", JHEP 07 
(2014) 079, doi : 10.1007/JHEP07 (2014) 079, arXiv: 1405.0301, 

[63] A. D. Martin, W. J. Stirling, R. S. Thorne, and G. Watt, "Parton distributions for the FHG", 
Eur.Phys.J. C 63(2009) 189, doi : 10.1140/ep jc/sl0052-009-1072-5, 

arXiv:0901.0002 

[64] M. R. Whalley, D. Bourilkov, and R. G. Group, "The Fes Houches accord PDFs (FHAPDF) 
and FHAGFUE", (2005). arXiv:hep-ph/0508110. 

[65] P. Artoisenet et al., "A framework for Higgs characterisation", JHEP 11 (2013) 043, 

doi:10.1007/JHEPll(2 013)043, arXiv:1306.6464 

[66] R. Frederix and S. Frixione, "Merging meets matching in MG@NFO", JHEP 12 (2012) 
061, doi: 10.1007/JHEP12 (2 012)0 61, arXiv: 120 9.6215 

[67] T. Sjbstrand, S. Mrenna, and P. Skands, "A brief introduction to PYTHIA 8.1", Comp. 

Phys. Comm. 178 (2008) 852, doi : 10.1016/j . cpc . 2008.01.03 6. 

arXiv:0710.3820 

[68] R. Gorke and T. Sjbstrand, "Interleaved Parton Showers and Tuning Prospects", JHEP 03 
(2011) 032, doi : 10.1007/JHEP03 (2011) 032, arXiv: 1011.175 9, 

[69] S. Bolognesi et al., "Spin and parity of a single-produced resonance at the FHG", Phys. 
Rev. D 86 (2012) 095031, doi : 10.1103/PhysRevD . 86.095031, arXiv: 1208.4018 

[70] T. Gorbett, O. J. P. Eboli, J. Gonzalez-Fraile, and M. G. Gonzalez-Garcia, "Determining 
Triple Gauge Boson Gouplings from Higgs Data", Phys. Rev. Lett. Ill (2013) 011801, 

doi:10.1103/PhysRevLett.111.011801, arXiv:1304.1151 

[71] G. Degrande et al., "UFO - The Universal FEYNRUFES Output", Comput. Phys. Commun. 
183 (2012) 1201, doi:10.1016/j. cpc .2012.01.022, arXiv :1108.2040 

[72] G. Altarelli, T. Sjbstrand, and F. Zwirner, "Physics at FEP2: Vol. 1", GERN Report 
GERN-96-01-V-1,1996. doi : 10.5170/CERN-1996-001-V-1 






25 


[73] CMS Collaboration, "Measurement of the W+W cross section in pp Collisions at = 7 
TeV and limits on anomalous WW 7 and WWZ couplings", Eur. Phys. J. C 73 (2013) 2610, 

doi:10.1140/epjc/sl0052-013-2610-8,arXiv:1306.1126 

[74] ATLAS Collaboration, "Measurement of WZ production in proton-proton collisions at 

= 7 TeV with the ATLAS detector", Eur. Phys. }. C 72 (2012) 2173, 
doi:10.1140/epjc/sl0052-012-2173-0,arXiv:1208.1390 

[75] ATLAS Collaboration, "Measurement of production in pp collisions at \/s = 7 

TeV with the ATLAS detector and limits on anomalous WWZ and WW 7 couplings", 
Phys. Rev. D 87 (2013) 112001, doi : 10.1103/PhysRevD. 87.112 0 01, 
arXiv:1210.2979, [Erratum: doi : 10.1103/PhysRevD .88.079906|. 

[76] CMS Collaboration, "Measurement of the W+W“ cross section in pp collisions at = 8 
TeV and limits on anomalous gauge couplings", (2015). arXiv:1507.03268 


A Appendix 





26 


A Appendix 


Table 1: Values of the pp —> H —> 77 differential cross sections as a function of kinematic 
observables as measured in data and as predicted in SM simulations. For each observable the 
fit to niryj is performed simultaneously in all the bins. Since the signal mass is profiled for each 
observable, the best fit oth varies from observable to observable. 
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